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1. Introduction 

Weak matrix elements of heavy-light mesons 
are of fundamental importance in particle physics 
since they enter determinations of some of the 
least known CKM matrix elements, and of the 
B — B and D — D mass splittings. The simplest 
matrix element to study is the leptonic decay con- 
stant of a pseudoscalar meson, fp. It is defined: 

< 0\A^\P >= fpPf,, 

where |P > is a pseudoscalar meson with 4- 
momentum p^, and is the axial current. 

The feasibility of using the lattice technique 
to calculate fp is now firmly established. The 
general trend in lattice calculations of fp is now 
towards a greater understanding of the system- 
atic effects entering the calculation. This has 
been made possible through smaller statistical 
errors uncovering systematics which were previ- 
ously hidden. Thus, the effects of, e.g., quench- 
ing, different choices of interpolating operators 
etc., can be studied. Systematic effects will be 
a focus throughout much of this review. 

The plan of this review is as follows. The next 
section overviews the three methods currently 
employed by the lattice community to calculate 
fp. Recent results from each of these methods 
are then presented in the subsequent sections. 
The present status of lattice calculations is sum- 
marised in the conclusion. 

Unless explicitly stated, all results are to be 
taken as "preliminary" . 
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2. Overview of Methods Used 

Heavy Quark Effective Theories (HQET) (see 
e.g. are invaluable tools in the study of the 
spectrum and decays of systems involving one or 
more heavy quarks. In this theory, the QCD ac- 
tion is systematically expanded in terms of the 
inverse heavy quark mass, mg. 

HQET is extremely powerful in lattice calcula- 
tions as well as in the continuum. It provides a 
means of overcoming the problem of simulating 
quarks with mass greater than the inverse lattice 
spacing. In general, a lattice calculation using 
HQET proceeds much as in conventional lattice 
calculations, but with all heavy quark propaga- 
tors calculated using HQET. Thus a heavy- light 
pseudoscalar meson correlation function in the 
conventional approach would be, e.g. 

C{t) = ^ < Ao(x)Ao(0) > 

= ^Tr{5,(x,0)75^Q(0,a;)75} 

X 

where Sq is the light quark propagator calculated 
using the traditional discretised Dirac equation, 
and the heavy quark propagator, Sq, is calculated 
using the discretised HQET action up to a chosen 
order in l/mg. 

In this review I will be discussing results ob- 
tained with Sq calculated using 

• the zeroth order HQET (also termed the 
'static' case) 

• HQET to some finite (typically the first) 
order (also termed NRQCD). 
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Figure 1. Representation of the range of validity 
of the three methods discussed in this review. 



I will also discuss results with Sq determined 
from 

• the conventional lattice Dirac equation 
(which I will term the "conventional" ap- 
proach). 

The next three sections discuss recent results 
from each of the above approaches in turn. The 
lattice actions or S'q's used in each case are de- 
fined in the corresponding sections. 

Some general comments regarding the applica- 
bility of the three approaches are helpful at this 
stage. One can represent the range of validity 
of each of the above methods graphically in fig. 1 . 
Here fpy/Mp is plotted against Mp. Note that in 
HQET it can be shown that, up to unimportant 
logarithms, 

jp^fu^ = constant + Oil/Mp). 

The static method is valid only for mg = oo and 
serves as a reference point for extrapolations in 
l/niQ. The NRQCD approach |,| is vahd for 
l/{niQa) « 1, and, at 0{l/mQ) its range of 
validity should include masses down to around 
Mb- The conventional approach, for present lat- 
tice parameter values, is believed to be accurate 



for meson masses of around Md and less. In 
fact, this statement assumes that some improve- 
ment scheme such as the SW action (hereafter 
termed the "Clover" action) and/or the "Fermi- 
lab" (also termed the "Heavy" ) formalism is 
implemented. If no such improvement is under- 
taken, then the conventional approach fails for 
Mp > Mo- 
Taken together, these three approaches span 
the entire range of mg , from infinity to the charm 
quark mass, nic and below. Clearly, a necessary 
condition for the success of lattice calculations 
is that these three methods agree in the shared 
regions, and that they provide a continuous func- 
tional behaviour for the quantities measured. It 
is one of the purposes of this review to study the 
status of this consistency check. 

To set the notation, lattice calculations of 
fp determine the dimensionless, unrenormalised 
quantity (p"^ corresponding to fp. It is defined 
(using an obvious notation): 

/p^Mp = Z^""(^#a-3/2. (1) 

(Note other notation for (p'^ appearing in the lit- 
erature includes V^Zl and v^/s-) I have chosen 
the superscript ^ to denote unrenormalised, di- 
mensionless lattice quantities. 

3. Static Results 

The static action has a very simple form and 
allows the heavy quark to propagate only in time, 
and not in position space. This action can be 
inverted to give the following form for Sq : 

Sq{x,0) = Vsit.) S{x) e{t,) 

where Vs{tx) is the Polyakov line from the point 
(f, 0) to {x,tx). 

3.1. Smearing 

The history of the static approach on the lattice 
has been plagued by problems associated with its 
poor signal to noise ratio. The reason for this 
effect is now well understood in terms of the vari- 
ance of the correlation function picture An 
approach to circumvent this problem is to intro- 
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duce smeared interpolating operators ||l2| . These 
can be expressed as 



{x, t)^Yl Q^y^ thol5q{z, t)iA 



y\ 



A)- 



Aq still couples to the ground state because it has 
the same quantum numbers as the local current 
Aq, but using it to define correlation functions 
can result in a smaller overlap with the excited 
states. 

There is a wide variety of forms of the smear- 
ing function ip presently in use. The most simple 
form of smearing functions is 'cubic' defined [p^ : 

for all <i3/2 ; i=l,2,3 
= otherwise. (2) 

These have the advantage that they are easiest 
to code, especially on parallel machines, but have 
the disadvantage that they require a gauge fixing 
procedure and are a rather brutal approximation 
to the physical wave function. This second point 



will be discussed more in sec 3.5 



Another popular choice of smearing functions 
is based on the non-relativistic quark model | p^ . 
Here, the heavy quark potential of the gauge con- 
figurations in the simulation is used to obtain the 
quantum mechanical wave function. Its advan- 
tage is that the physics is being used to dictate 
the smearing. 

The Wuppertal collaboration have proposed 
various smearings . The advantage in these 

cases is that the smearing is relatively quick, and 
is gauge invariant. 

My favourite smearing is "MOST" (Maximal 
Operator Smearing Technique) [0 . In this tech- 
nique the set of operators is defined where the Q 
and q fields have every possible separation (which 
is independent under cubic symmetry) . The two- 
point correlation between all members of this set 
is measured, and an simple analysis [p"8|jig| ] can be 
used to extract the ground state properties. The 
advantage of this method is obvious: because it 
uses all possible relative separations as a basis, all 
other smearing methods are a subset of this ap- 
proach. Its use in light-light physics should there- 



fore be investigated. The disadvantage is that it 
becomes memory intensive for large lattices. 

Note that in order to extract the local matrix 
element, /p, local-smeared, as well as smeared- 
smeared correlation functions need to be mea- 
sured. In practice it is best to place the local 
current at the sink rather than the source, since 
a better signal to noise ratio results ||2^,|l^ . 

3.2. Z^^" 

In order to define the physical fp a choice needs 
to be made for the definition of the renormal- 
isation constant Z^*^" see eq.(^. Pcrturbative 
calculations j2^ show, in the case of the static- 
Wilson current, 

ZftTAa^) - (l-.18992)(l+.0127.g2Zo5(a2rn|)).(3) 

For the static-Clover case, 0.189 is replaced by 
0.144 Ig^. Note that the first factor is essen- 
tially the lattice to continuum matching within 
the effective (static) theory, and the second is the 
matching between the effective and full theory in 
the continuum. 

Using the Kronfeld-Lepage-Mackenzie formal- 
ism [^,0 leads to the following definition of 
for the Wilson action. 



en 
t 



7Ren 



X (1 -0.137g2) 



X (1 + ^m2TgHog{a^m\y). (4) 

Here, the renormalised couplings are = &/ fi/u^ 
and Kq = UQKq. uq will be defined later. The 
first factor in eq.(Q) is the rescaling of the light 
quark field from the traditional, Wilson y^2Kq to 
the KLM prescription. There is however an am- 
biguity in the above definition of Z^'^^. This is 
because one is free to define uq as, e.g. the aver- 
age plaquette, < Upiaq >^^'^^ or l/8/«c (where Kc 
is the critical value of the hopping parameter). 

Table |l| displays the values of Zf^'^^ obtained 
with the above definitions of Zffjj:. As can be 
seen from the table the difference between the 
various determinations of Z^"^"^ is 0(30%) and it 
is 0(5%) even within the KLM formalism. Obvi- 
ously, this means that the overall systematic er- 
ror in Z^*"* has a contribution of this order from 
the uncertainties in Z^*^". This error from the 
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Table 1 

ZjJ^" obtained using different definitions. Col- 
umn 2 uses the "naive" definition of = 6//3 (see 
eq.(^). Columns 3 and 4 use the KLM formula- 
tion (at the chiral limit) with the uq definitions 
given by 1/(8kc) and < Upiaq >^^^ respectively 
(see eq.(ll)). 



yHen 
^stat 



6.0 
6.1 
6.2 
6.3 
6.4 



0.83 
0.83 
0.83 
0.83 
0.83 



0.62 
0.64 
0.66 
0.67 
0.68 



0.59 
0.61 
0.63 
0.64 
0.65 



determination of Z^f^ has now become one of 
the dominant errors in the lattice calculations of 
/g"*. The way around this problem is to use 
a non-pcrturbative definition of Z^"^" which are 
becoming available which circumvent the above 
problems completely [Q. 

3.3. Continuum Limit 

The possible variation of fp in the static limit 
with the lattice spacing, a, has been long dis- 
cussed (see, for example, |2^,|l^,|lJ] This so-called 
non-scaling behaviour has been proposed as a 
mechanism for reducing the lattice prediction for 
fp by as much as 30%. In principle this a de- 
pendence is simple to study. One simply calcu- 
lates fp on the lattice using a number of /3 values 
and extrapolates the results to a = 0. In prac- 
tice this method is difhcult to implement due to 
the increasing statistical errors of the data points 
closest to the point of interest, a — Q. 

Recently a method has been proposed to anal- 
yse the scaling behaviour of fp, and indeed any 
other dimensionful lattice quantity p6| . Instead 
of studying the scaling of the final quantity, /p, 
the scaling of each of the three factors in eq.(|l|) 
which define fp: Z^^^; and are deter- 
mined separately. Each of these three factors are 
functions of g^. The question to be answered is: 
does the g^ dependence of these factors cancel in 
the overall product? The analysis performed in 
Pql suggests that, within present statistical er- 



rors, fp does scale but only for /? ^ 6.0. This 
means that analyses which use a linear fit in a 
and data with (3 ^ 6.0 are biasing their contin- 
uum extrapolation to a smaller value. This is due 
to the levering effect of the non-scaling, and typi- 
cally higher data for /3 ^ 6.0, which, furthermore, 
typically have smaller statistical errors. 

The analysis in ||2^ also showed that the dom- 
inant non-scaling effect was not due to higher 
orders in perturbation theory in g'^ = 6//3, but 
rather to 0{a) effects. A common 20% to 30% 
discretisation effect was found for the quantities 
studied, such as (p"^, Mp and f^^. 

3.4. Recent Results 

Recent results for /^*°* are shown in table ^ 
together with their references in the second row. 
A striking fact about these latest simulations is 
the appearance of simulations with dynamical 
fermions 28|. A discussion of the effects of 
this on fp appears in sees. 4.2,5.1 and in the 
conclusion. 

A brief description of each of the groups' work, 
together with a longer discussion on the results 
from the APE collaboration follows. I have cho- 
sen to present a relatively detailed discussion of 
the APE collaboration's work since it will not be 
presented elsewhere in the proceedings. 

The UKQCD Collaboration Q presented re- 
sults for both /f and Bg''* (see also They 
noted that the systematic error associated with 
^Ren jg j^Q^ ^ major uncertainty (see sec. p.2|) . 
A summary of their results is shown in table 
(The first error is statistical and the second sys- 
tematic.) 

The SCO collaboration ||l|,|3 used both the 
static and NRQCD implementation of the heavy 
quark propagator, and their results will be dis- 
cussed in detail in the NRQCD section ^. 

The MILC collaboration |Q combine an anal- 
ysis of both the static and conventional ap- 
proaches, and a discussion of their work will ap- 
pear in the latter section (sec. 5.1). 



3.5. APE Results 

The APE collaboration has recently carried 
out a high statistics calculation of Z^*"* in the 
quenched approximation p2] . This work is a con- 
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Table 2 



Summary of recent results using the static approach. 'Q' signifies Quenched simulation, 'D' Dynamical. 
See text for detailed comments. All numbers should be considered preliminary with the exception of [30]. 



Collaboration 

Reference 


UKQCD 
p9|,p0| 


SCO 


MILC 
rFToli 

& 


APE 
|j3z| & see text 


Lattice Parameters 













& 2 


& 2 





P 


6.2 


5.6 D 


5.7^ 6.5 Q 


6.0,6.1,6.2,6.4 






6.0 Q 


5.4-> 5.7D 




Volume 


243 X 48 




< 323 X 100 


< 243 X 64 




60 


100 D 


O(IOO) 


< 400 






35 Q 








Clover 


Wilson & Clover 


Wilson 


Wilson & Clover 




Static 


Static 


Static 


Static 




& NRQCD 


& Conventional 




Results 

[MeV] 


Of; (^+18+28 
^OO_20_27 


see Section^ 


see Section || 


240(30) 



/|*;V/f"* i.iet^ " " 1.17(4) 

Mbs - Mb [MeV] 87j:i^t^2 ^ " 81(10) 



tinuation of a long program of work in this area 
ISS^pHl . As discussed in sec. 3.1 , the cubic smear- 



ing function was used. The methods used in the 
past for extracting the local matrix element (p"^ 
from the smeared correlation functions are dis- 
cussed in [Q. One problem with these methods 
is that they require the selection of the "best" 
cube size Lg (see eq.(j^)). A number of techniques 
were developed to make this choice as free from 
bias as possible [^,^. However, the problem re- 
mained that cubes are available at only discrete 
sizes, and that as the light quark mass was var- 
ied, the meson's variation in size could not be 
adequately tracked by the cubes. For example, 
at (3 = 6.4 the best cube size for the heaviest 
of the light quarks studied was clearly Lg = 15, 
but for the lightest quark was somewhere between 
Lg = 15 and Lg = 17. 

A solution to this problem has already been 
proposed by a number of authors (see e.g. [p^ ): 
one includes various cube sizes altogether in the 
one fit so that the fit is not constrained by a 
single cube size. Thus, in the example considered 
above for (3 — 6.4, the new enlarged is 



E E 



- E 

Ls = 15,17 t i=SS,LS 



In this definition, both the Lg = 15 and Lg — 17 
cubes are included in the as well as the usual 
sum over the times t in the fitting window. The 
sum i — SS, LS includes the contributions from 
the smeared-smeared and the smeared(at source)- 
local(at sink) correlation functions. The fitting 
function /^^'*(t) is chosen appropriately; for ex- 
ample. 



/ 



Ls = 15,i=LS 



it) 



15) 



e-*^* X 



{l+v{LMl5)e-^') (5) 



where a two state fit is performed and the ?7's 
and A are the overlap and mass gap of the first 
excited state respectively. 

The APE Collaboration has re-done its analysis 
for its whole data set at P — 6.0,6.1,6.2 and 6.4 
using the above method. As an example, fig. 2 
shows the results of the new method together with 
the previous analysis in the case of /3 = 6.4. As 
can be seen, the use of the new method results 
in an increased slope in the chiral extrapolation. 
Thus the chiral ratios and differences, /sgV/s"* 
and Mbs — Mb will increase. The overall effect 
is not large, but, because the statistical errors are 
getting smaller, systematic effects such as this are 
becoming important. 
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Figure 2. APE results for (fy^ as a function of 
the quark mass I/kq — 1/kc- The data shown is 
for /3 = 6.4 with the Clover action for the light 
quarks. The plot shows the results of the old 
method for = 15 ("+") and 17 ("□"), and the 
new method ("O") (see text). Points from the 
new method have been shifted horizontally to aid 
clarity. 



In fig. 3, the values for /^*°* at various j3 values 
are plotted against the lattice spacing a using the 
new method. In this plot, the rho mass was used 
to set the scale. A preliminary fit to a constant 
(see sec.O) results in /f = 240(30) MeV. Val- 
ues of /^//f^* = 1.17(4), and Mbs - Mb ^ 
81(10) MeV (compared to 96(6) Mev from ex- 
periment) are also predicted. All errors quoted 
here include both the statistical and a — > un- 
certainties. Finite volume effects are believed to 
be negligible for these lattices ||l^ . 

4. NRQCD Results 

The NRQCD action used in current simulations 
is correct to 0(1/toq) 

C = Q^ (-A + ^ + tt^'^-b] Q. 

[ 2mQ 2mQ J 

This action can be simply inverted to obtain the 
heavy quark propagator. In addition to the 1/toq 
terms which appear in the action, there are 1/toq 
terms which appear in the definition of the cur- 
rents p6t . These must be included in order to 
make the calculation correct at this order. All 



Figure 3. APE results for /^*°* for the various 
simulations as a function of the lattice spacing a 
(set from the rho mass). From right to left the 
data is for Clover /3 = 6.0, Wilson /3 = 6.1, Clover 
(3 = 6.2, Wilson P = 6.4 and Clover /? = 6.4. 
Horizontal errors are suppressed for clarity. 



the NRQCD work presented at this conference 
took these terms into account. 

The NRQCD approach is valid for heavy-light 
meson masses down to around the B-meson. Un- 
fortunately it is not valid at the charm mass 
where higher order terms in the l/mga expan- 
sion become relevant ||3^ . In theory, moving to a 
larger a value (smaller /3) would mean that these 
terms become less relevant. In practice however, 
this would only increase the 0{mqa) effects to an 
extent where the results would be unreliable in 
any case. 

Table I summarises the data presented in this 
conference with the references appearing in the 
second row. The following sections detail the re- 
suhs of the SCO (SCRI- Glasgow-Ohio) collabo- 
ration. In the first, the quenched approximation 
was used, and in the second dynamical quarks 
were included. See |Q for details of the Ken- 
tucky calculation. 

4.1. SGO Results: Quenched 

The SGO collaboration's quenched data ||3l|| 
are run using the Clover action for the light 
quarks: both the tadpole unimproved and im- 
proved (i.e. c — 1 and ~ 1.4 respectively, where 
c is the coefficient of the Clover term in the ac- 
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Table 3 

Summary of recent results using the NRQCD approach, 
should be considered preliminary. 



See text for detailed comments. All numbers 



Collaboration SGO 

Reference [Bll 



SGO 



Kentucky 
^381 



Lattice Parameters 



P 

Volume 

Ncfgs 
Sq 





6.0 

16^ X 48 
- 50 
Clover 
NRQCD to 0{l/mQ) 
& Static 



2 

5.6 

16^ X 32 
100 

Wilson & Clover 
NRQCD to Oil/mg) 
& Static 





6.0 

20^ X 30 

32 
Wilson 
NRQCD to Oil/mg) 



Results 

fs [MeV] 

£stat I fstat 
JBs IJB 



160(40) 



- 190 

1.18(3) Clover 
1.28(4) Wilson 



1.15: 



tion). Two values of Kq were used which strad- 
dle the mass of the strange quark. Four values 
of the NRQCD quark mass were used. The re- 
sults for (jy^ at the chiral limit versus l/Mf are 
plotted in fig. 4 for the c — I case together with 
the static point. Fitting only the NRQCD points 
shows that they can be made to extrapolate to the 
static point only if a quadratic term, C2/Mp, is 
included in the fit. If a linear term only is allowed, 
then the static point seems too high. However, it 
may be that with the inclusion of more mg values 
closer to the static point , th is apparent discrep- 
ancy disappears (see sec.O). Also, the c = 1.4 
data (not shown) has a better agreement between 
the NRQCD data and the static point ]3l|] . 

A preliminary calculation of the renormalisa- 
tion constant relevant for the NRQCD-light axial 
current has been performed p9|| . Using this value, 
the SGO collaboration obtains fp at their 4 val- 
ues of mg. From these, I obtain the preliminary 
estimate of Jb which appears in table ||. 

4.2. SGO Results: Dynamical 

Dynamical simulations involving NRQCD were 
presented by the SGO Collaboration in last year's 
conference |Q. This year's results are en- 
hanced by the inclusion of more values of mg, 
by a further analysis of the ©(l/mg) terms that 
contribute to fp, and by the use of the tad- 
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1/Af 



# 



Figure 4. Quenched SGO results for cj)'^ as a 
function of the inverse heavy-light meson mass 
l/A/p# at the chiral limit [31]. The results of a 
linear and quadratic fit in 1/Mp = \/{Mpa) to 
the NRQCD data points is also shown. 



pole improved Clover action for the light quarks. 
Again the preliminary values for the renormali- 
sation constant were used For the Wilson 

data only two values of the light quark mass were 
available meaning that chiral extrapolations were 
not robust. The Clover data has three light quark 
values. 

In fig. 5 a plot of Z^^'cj}* = Jp^/M^a^/'^ (see 
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eq.(Q)) versus l/Alf is displayed. In this graph 
both Wilson (at Kg — 0.1585) and Clover (k^ = 
0.1385) data are plotted. The Kg values are cho- 
sen so that the corresponding pion masses agree. 
Note that in order to make the comparison, the 
appropriate renormalisation factor for both ac- 
tions has been included. The curves are cubic fits 
in l/M-jt to the NRQCD data only. The static 
points are also plotted, but not included in the 
fits. 

There are three important points to be made: 

• The NRQCD data points smoothly extrap- 
olate to the static point. Thus the inclusion 
of more data points with larger mq seems 
to have resolved the discrepancy that ap- 
peared earlier po| . 

• The slope of 4> against 1/Mp is appar- 
ently increasing as Mp — > oo. This signi- 
fies the presence of terms of order l/mg. 
However, the calculation is correct only to 
0{l/mQ). Therefore, strictly speaking, any 
non-constant behaviour of this slope at fi- 
nite niQ is not a real prediction of these 
calculations. In practice though, the be- 
haviour seen in fig 5 matches well with data 
from conventional simulations at around 
the D-meson (see sec. |^). So presumably it 
true is that the present simulation does not 
suffer from neglecting terms of order l/rriQ. 
Note also that from the figure there is a 30% 
l/Mp correction to /b compared with the 
static value. 

• The Wilson data is statistically lower than 
the Clover data. The obvious explanation 
is that this is a symptom of 0(a) effects 
spoiling the Wilson data (These sys- 
tematics could enter either in (j)'^ or in a~^.) 
It could also be that the one-loop pertur- 
bative calculation of Z^"^" is inadequate, 
and that a non-perturbatively defined Z^'^"' 
would remove the disagreement. 

In the SCO analysis the contributions of the 
three ©(l/mg) terms to cj) were extracted. Ex- 
pressing these as Ci/Mp they found that 

|cli2| > \Ca,D\ » \Ca.B\ 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 



1/M* 

Figure 5. Dynamical SCO results for Z^'^^cj)'^ as 
a function of the inverse heavy-light meson mass 
1/M* [27]. The upper set of data uses the Clover 
action for the light quark (with Kg = 0.1385), 
and the lower set uses the Wilson action (with 
Kg = 0.1585). The results of cubic fits in 1/Mp 
to the NRQCD data points are also shown. 



That the hyperfine component, Co-.s is smallest 
fits our expectations since this term breaks both 
the spin and flavour symmetry present in the 
niQ = oo hmit ||l|,|3- 

Performing a chiral extrapolation and using a 
nominal value of a^^ = 2 GcV, one obtains the 
values for fp and fp/fBs which appear in ta- 
ble || H]. This (rough) fp value IS m nice 
agreement with that obtained by the quenched 
NRQCD analysis of js^ implying that the effects 
of dynamical quarks are at, or below, the level of 
statistics. The difference between the two actions 
for the chiral ratio is unexpected since systematic 
effects such as 0{a) are expected to cancel in this 
ratio. 

5. Conventional Results 

The "conventional" formalism of the heavy- 
quark propagator is simply that obtained from 
the lattice version of the Dirac equation, i.e. the 
same lattice action as for the light degrees of free- 
dom. However, potential problems exist when 
iriQa <; 1. There are two approaches to reduce 
discretisation effects. The first is to use an im- 
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proved action, such as the Clover action W, which 
has all terms 0{mQa) explicity removed [|42| . The 
second is to use the "Fermilab" formulation of |^ 
||. These two choices are not mutually exclusive, 
and while all groups use the second technique, 
one group |Q uses both. 

Table | summarise the data to appear recently 
with the references in the second row. In the fol- 
lowing subsections a brief review of each of the 
four groups' work is presented. 

5.1. MILC Results 

The MILC collaboration's results Q distin- 
guish themselves because they are the only simu- 
lations in the conventional approach that (i) have 
dynamical fermion runs, (ii) have a static point 
to help constrain their fits in 1/Mp, and (iii) use 
the hopping parameter expansion to give them a 
wide range of kq values at very little cost [Q. 
They have also studied finite volume effects, and 
performed a continuum extrapolation. (This last 
point will be discussed more later.) 

In fig. 6 the quantity (j) is plotted against l/Mp, 
both in physical units. The curve shown is a fit 
to the points marked with a cross (i.e. in this 
case, points with Mp < Mp,)- By varying the 
fitting window in Mp, an estimate is obtained of 
the error due to the mqa effects which escape 
the Fermilab redefinitions. Another estimate of 
this error would be to fit the standard Wilson 
data (without the Fermilab redefinitions) for, say, 
Mp < Md, and compare it with the results from 
Fermilab data for the same Mp range. 

Fig. 7 shows the continuum extrapolation of fp 
with used to set the scale. Both quenched 
and dynamical data are shown. The line shown 
is a linear fit to all the quenched data (i.e. 5.7 < 
/3 < 6.5). It is probable that this procedure suf- 



fers from the same problem discussed in sec. 3.3 
where static data for (3 ^ 6.0 was found not to 
scale, within statistical errors. If this scenario 
is correct, then including the data with /3 ^ 6.0 
in a linear fit skews the continuum extrapolation 
downwards. Fig. 7 seems to confirm this assertion, 
though the effect is small. This would explain at 
least part of the difference between the MILC re- 
sults and the other data shown in table ^. 

The dynamical data points in the figure show a 



/9 = 6.3, Lighter heavies in fit 



Ki used in fit 




0.0 0.5 1.0 

1 /Mp (GeV^^) 

Figure 6. MILC results for (/) as a function of 
1/Mp [28]. The data points used in this particu- 
lar fit are denoted with a cross (see text). 



fg; lighter heavies; scale 



quenched 



> 



4+ 




150 



a (GeV)"^ 

Figure 7. MILC results for fp as a function of 
lattice spacing a [28]. The fit is to the quenched 
points only (see text). 



clear discrepancy at around a — 0.5GeV~^ com- 
pared with the quenched data. This can be used 
as an estimate of the quenching errors which ap- 
pears to be considerable at this stage [|2^ . 

The results of the MILC collaboration's work 
is displayed in table IJ. The three errors stated 
are: (i) statistical, (ii) various systematic, and 
(iii) quenching. The systematic errors are a com- 
bination of the uncertainties in the various fits, 
the mqa effects, finite volume and finite a errors. 



10 



Table 4 



Summary of recent results using the conventional approach. 'Q' signifies Quenched simulation, 'D' 
Dynamical. See text for detailed comments, and a description of the various errors. All numbers should 
be considered preliminary. 



Collaboration 


MILC 


LANL - Ohio 


JLQCD 


FNAL 


Reference 


PI 


Washington 

|45 


0] 


m 


Lattice Parameters 










& 2 













5.7^ 6.5 Q 
5.4^ 5.7 D 


6.0 


6.1,6.3 


5.9 


Volume 


< 323 X 100 


323 X 64 


< 323 X 8Q 


163 X 32 


^cfgs 


- 100 


170 


~ 100 


100 


Sq 

Sq 


Wilson 


Wilson 


Wilson 


Clover 


Fermilab- Wilson 


Fermilab- Wilson 


Fermilab- Wilson 


Fermilab-Clover 


& Static 








Results 










Id [MeV] 
fB [MeV] 

fBs 


182±3(9)(22) 
198±5(10)(19) 
151±5(16)(26) 
169±7(14)(29) 


229 ±7(-)(+12)(li«)(7) 
260±4(-5)(+15)(±^4)(8) 


214 ± 10(25) 
248 ± 7(29) 
221 ±15(26) 
244 ±8(28) 


220^5 
(23913) 
188t^) 
(207^3) 



5.2. LANL-Ohio- Washington Results 

The LANL-Ohio-Washington col- 

laboration have new results using the Fermilab - 
Wilson formulation ^] . Since they have no static 
point, and are working at a moderate P value, no 
estimate of Jb was given. Of particular interest 
is their comprehensive study of the systematics 
involved in setting the strange and charm quark 
masses, and in the definition of Z^'^". This is 
outlined below. 

• The hopping parameter corresponding to 
the strange quark mass, Ks, is defined in 
three different ways: by fixing Mj^/M^; 
M'^/Mp-, and AI^/Mp to be equal to their 
physical values. These estimates corre- 
spond to a surprisingly large ~ 20% spread 
in the values of the strange quark mass. 

• The charm quark is set by fixing the mass of 
a meson containing a c-quark to its experi- 
mental value. The lattice estimate of these 
heavy-light meson masses is defined using 
either the pole mass (obtained from the 
exponential decay of the correlation func- 
tions), or by the lattice dispersion relation. 



• A careful study of the systematic effects 
entering the definition of Z^'^" was under- 
taken. Here different values of the Lepage- 
Mackenzie q* were used [Q, and two defi- 
nitions {uo = 1/8kc and < Upiaq >^^^) were 
used to redefine the coupling. 

The spread of decay constant values using each 
of these methods gives estimates of the corre- 
sponding systematic effects. The final results for 
and fu^ are shown in table ^. The five errors 
listed are due to: (i) statistics, (ii) strange quark 
determination, (iii) charm quark determination, 
(iv) Z^'^", (v) setting the scale. The importance 
of this work is that it shows that each of the three 
errors outlined above are at least of the same or- 
der as the statistical and scale errors (which are 
normally assumed to be the dominant errors). 

5.3. JLQCD Results 

A status report of the ongoing analysis by the 
JLQCD collaboration was presented at this con- 
ference 0] . They have data at two {3 values and 
use a quadratic fit in 1/Mp to determine their 
values for /s, fo etc. Their results, taken from 
the /3 = 6.3 dataset are given in table |. The p 
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mass was used to set the scale. The first quoted 
error is statistical, and the second is due to the 
scale. More accurate results will be obtained from 
this collaboration in the near future. 

5.4. FNAL Results 

The FNAL group are simulating the Clover ac- 
tion in the Fermilab formalism at (3 = 5.9 p3[ . 
They have preliminary results which are shown 
in table ^. (Statistical errors only are shown.) 
The strong point of this work is that it uses the 
Clover action (i.e. complete with the rotations 
of the quark fields) and that therefore the results 
should be free of 0{mQa) effects. Values for the 
decay constants in the table use to set the 
scale. stnd Jbs are quoted in parentheses, 

since the strange quark value used in their deter- 
mination is nominal. The FNAL group plan to 
simulate at other /? values to check the scaling of 
these quantities. 

6. Conclusions 

One of the main themes of this review is that 
systematic effects now clearly dominate statistical 
errors in lattice calculations of fp. In many ways 
this is a very desirable situation since it allows 
their careful callibration. I list the main system- 
atic uncertainties and their possible solutions. 

• Z^'^'^ varies by as much as '-^ 30% (see sees. 
3.2 and 5.2). This variation can be sig- 
nificantly improved by using the clover ac- 
tion. The non-perturbative value of Z^*^" 
will clear up this uncertainty [ p4| . 

• The continuum limit of fp can be studied 
by an a ^ extrapolation. In the static 
case, only data with /? ^ 6.0 should be in- 
cluded in this fi t (fo r present levels of statis- 
tics, see sec. ^.3|) , otherwise the a — 
value will be biased downwards. Presum- 
ably a remnant of this effect survives in the 
NRQCD and conventional cases. Clearly 
0{a) effects should be reduced by the use 
of the clover action. 

• Dynamical {np — 2) results seem to pro- 
duce higher values for fp for both the 
NRQCD (sec. and conventional cases 



(sec. 5J_). Whether this effect is invariant 
under further investigation will be a main 
source of research over the next couple of 
years. A more subtle question is the extent 
to which np = 2 data truly refiects the real 
npKi3 world @,|§. 

• The choice of physical quantity to set a im- 
pacts upon the final value of fp. Again, 
short-term, this should be fed into the sys- 
tematic uncertainties, but longer-term this 
should be understood in terms of quenching 
and/or 0(a) errors. 

• The chiral extrapolation of methods which 
use a fixed smearing for all light quark val- 
ues should be questioned (see This is 
not a large effect for fp but can become sig- 
nificant for chiral ratios/differences such as 



fBs/fs and Mbs 



M 



B- 



• Setting the hopping parameters correspond- 
ing to the physical strange and charm 
quarks also introduces an error which 
should be included (see sec. 5.2). 



The wide range of Mp values presently cov- 
ered by the three methods allow an important 
consistency test to be made: Do the different ap- 
proaches agree with eachother? Studying the fs 
values in tables ^,|| and |4[ together with the ver- 
sus 1/Mp plots suggests that they do. Note that 
there is still a fair degree of 'slop' in the data 
due to the fact that different groups have used 
different definitions of Z^'^^'' and chosen different 
quantities to set the scale. Once this is tightened 
up, this consistency test can be made more defini- 
tive (e.g. by using all three methods to determine 
fp on the same configurations). 

Using the data in the tables I estimate the 
following "global" lattice averages: fo ~ fs = 
200MeV ± 20%. The error bar includes my esti- 
mate of all uncertainties. 
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